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I— I 10 years from now medium baseline reactor experiments will attempt to determine the neu- 
p I trino mass hierarchy from the differences {RL + PV) between the extrema of the Fourier 
transformed neutrino spectra. Recently Qian et al. have claimed that this goal may be 
fH impeded by the strong dependence of the difference parameter RL + PV on the reactor neu- 
trino flux and on slight variations of |AM|2|. We demonstrate that this effect results from 
a spurious dependence of the difference parameter on the very high energy (8+ MeV) tail 
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O of the reactor neutrino spectrum. This dependence is spurious because the high energy tail 

00 depends upon decays of exotic isotopes and is insensitive to the mass hierarchy. An energy- 

^ dependent weight in the Fourier transform not only eliminates this spurious dependence but 

O in fact increases the chance of correctly determining the hierarchy. 
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This year the Daya Bay [1, 2] and RENO [3] experiments have demonstrated beyond any 
reasonable doubt that Oi^ is as much as an order of magnitude larger than had been suspected 
several years ago. This large value of ^13 implies that 1-3 neutrino oscillations may be 
observed at medium baselines, which we define to be 40-80 km. The medium baseline 
neutrino spectrum may then be used to determine the neutrino mass hierarchy [4]. Such 
experiments are now not only practical but indeed they will be performed within the next 
decade [5, 6, 7]. 

However at these baselines, due to a degeneracy in the high energy neutrino spectrum 
[8, 9, 10], a determination of the hierarchy requires a measurement of 1-3 oscillations at low 
neutrino energies E. As a result of the finite energy resolution of the detector and various 
interference effects [10] these low energy peaks are difficult to identify individually at medium 
baselines L. Nonetheless if the nonlinear energy response of a detector is well understood 
then one may measure the sum of the peaks by studying the k ~ |AM|^|/2 region of the 
L/E-Fouher transform of the neutrino spectrum [11]. The most popular variables for such 
a determination are the fractional difference RL between the deepest minima of the Fourier 
cosine transform and the difference PV between the deepest minimum and the highest peak 
of the Fourier sine transform [12, 13]. Although these two variables are somewhat degenerate, 
an improvement may be obtained by considering their sum RL + PV. 

A serious obstruction to this analysis, and thus to plans to measure the neutrino mass 
hierarchy at medium baselines, has been described in Ref. [9]. The authors observed that 
the combination RL + PV is very sensitive to the choice of model of the reactor neutrino 
flux and to variations of IAMI2I which are smaller than the precision to which this 

mass difference has been determined by MINOS [14]. While the observed shift appears to 
depend upon both IAMI2I and the hierarchy, from Fig. 4 of Ref. [9] it can be seen that the 
shift depends only upon the effective mass difference [15, 10] 

AM^s = cos\9u)\AMi^\ + sm\9i2)\AMl^\. (1) 

The neutrino flux from reactors is known poorly. The theoretical normalization has recently 
increased by about 3% [16] and the 6+ MeV flux has increased by an additional 3% [17]. 
The flux beyond about 8 MeV is not known at all due to its strong dependence upon decays 
of exotic isotopes [16]. Even worse, all of these theoretical fluxes are about 6% above the 
observed fluxes at very short [18] and 1 km [19] baselines. Thus the large sensitivity of 
RL + PV upon the poorly known fluxes and AM^g appreciably reduces the probability that 
a medium baseline reactor experiment can correctly determine the neutrino mass hierarchy. 

We will now explain the cause of this strong dependence. As the dependence of RL, 
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PV and RL + PV upon these parameters is virtually indistinguishable, for brevity we will 
consider only 

which is the fractional difference between two minima R and L of the Fourier cosine transform 
of the neutrino spectrum 

where the tree level neutrino inverse /3 decay cross section is [20] 
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a(E) = 0.0952 x 10"^^ cm^ ^ % — E^ = E - mn + (4) 
and the electron neutrino survival probability is 

Pee = Sin^(^i3) + C0S4(^12)C0S^(^13) + sin^(^i2)cOS^(^i3) + ^(Pl2 + ^13 + ^^23) 

P12 = sin2(2^i2)cos^(^i3)cos ( , Prs = cos2(^^i2)sin2(2^i3)cos ^'^^'J^ 
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Following Ref. [12], a 3%/\/E energy resolution is included by convoluting the observed 
energy spectrum with 

(E - 



""""^V 0.0018(Ee + me)MevJ • 
We use the neutrino mass matrix parameters of Ref. [5] . 

As was demonstrated in Ref. [10], the minima whose difference defines RL lie just on 
either side of A; = |AM|i|/2. These minima arise from the Fourier transform of P13 which 
is independent of the hierarchy, but the contribution of P23 provides a perturbation which 
makes the right (left) minimum deeper for the normal (inverted) hierarchy. 

The problem observed in Ref. [9] is that, depending upon the reactor flux model used, 
the transform of the unoscillated reactor flux ^{E)a{E)E'^ / itself may contribute peaks 
near k = |AM|^|/2 which interfere with those of P13 + P23 and so affect RL. While the 
cosine transform of the unoscillated flux ^{E)a{E)E'^ / L^ is independent of the neutrino 
mass splittings, the locations of the peaks of P13 + P23 are proportional to AM^g. This 
means that the relative phase between the Fourier transform of the unoscillated spectrum 
and that of P13 + P23 depends on the precise value of AM^^. As a result the oscillations in 
the Fourier transform of ^{E)a{E) lead to an AM^g-dependence in the quantity RL just 
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Figure 1: The cosine transforms of the unoscillated flux (black dashed curve) and the full 
Pi3 + P23 oscillated flux (red solid curve) are shown. As AM^^ varies, the P13 + P23 peaks 
move and so the interference between the two contributions to the cosine transform of the 
neutrino spectrum varies. The reactor fluxes used are those of the 1980's. 

of the kind observed in Ref. [9] using old reactor flux models. In fact, using the ^ss-jj g^^^ 
from Ref. [21], the ^^^Pu and ^^^Pu fluxes from [22] and the Gaussian approximated ^38^ g^^^ 
from Ref. [23] with the isotope ratios of Ref. [12] we find an oscillation in the unoscillated 
spectrum term in Eq. (3). Using this old model of the reactor flux, in Fig. 1 we compare 
the Fourier transform of the unoscillated term with that of the P13 + P23 term, which is 
sensitive to the hierarchy. One can see that the unoscillated term is periodic with the same 
wavelength as was observed in Fig. 4 of Ref. [9] , and thus the interference between these two 
terms oscillates as AM^^ varies, shifting the P13 + P23 peaks and so reproducing the effect 
reported in that note. 

Ref. [9] concludes that this strong dependence of RL upon the reactor flux means that a 
precise knowledge of this flux is desirable to determine the neutrino mass hierarchy at a short 
baseline experiment. Our conclusions differ. This difference arises from the observation that, 
as can be seen in Fig. 2, near k = \AM^-^\/2 the oscillations of the cosine transform of the 
reactor spectrum are governed by the highest energy neutrinos (8+ MeV) which, according 
to Refs. [8, 10], are not sensitive to the hierarchy. On the other hand in Fig. 2 we see that the 
effect is present in the case of both old and new fluxes if the spectrum is cut off at 8.5 MeV 
while it is negligible if the spectrum is cut off at 12.8 MeV. However since the quadratic and 
quintic fits to these fluxes are only reliable below 6.5 MeV and marginally reliable below 8 
MeV, there is no reason to trust a naive extrapolation to 12.8 MeV. 

As RL depends strongly on the spectrum between 8.5 and 12.8 MeV, which in turn is 
independent of the hierarchy, this high energy spectrum provides a nuisance parameter for the 
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Figure 2: The cosine transforms of the unoscillated flux is shown for numerically interpolated 
fluxes from the 1980's [21, 22, 23] (black dotted curve), for a quadratic fit to fiuxes from the 
1980's [23] and for quintic fits of the new fiuxes of Ref. [17]. The latter two are shown with 
cutoffs of 8.5 MeV (dashed curves) and 12.8 MeV (solid curves). The blue dashed curve 
corresponds to the quadratic fit fiux. The red and green solid curves, corresponding to 12.8 
MeV cutoffs, are close to zero. Therefore the interference effect is present if the cutoff is at 
8.5 MeV and but not if the fits are naively extrapolated to 12.8 MeV. This demonstrates 
that RL is sensitive to the neutrino spectrum above 8.5 MeV. 



determination of the hierarchy using RL. The solution suggested in Ref. [9] is to determine 
the spectrum precisely, however so few neutrinos are observed in this range that such a 
determination would be difficult, indeed the spectrum is not understood at the required 
precision even at the energies with high fiuxes [18]. Even if such a measurement were 
possible, then RL would still be likely to depend upon AM^g with a higher sensitivity than 
the mass determination at MINOS and probably at NOz^A, making a determination of the 
hierarchy at a medium baseline more challenging. 

Our solution is to replace RL and PV with indicators that are insensitive to the high 
energy neutrino spectrum, by providing an energy-dependent weight w{E) on the neutrino 
spectrum in the Fourier transform. As we saw in Fig. 2, a simple cutoff in the Fourier 
transform will amplify the spurious dependence. The weight needs to cut off the high energies 
gradually, with derivative scales much longer than |AM|]^|, so as to not itself introduce 
spurious peaks in the critical part of the Fourier transforms. One such choice of weight 
which we have found works quite well is a Gaussian 

^^^^ f J L\ o.osE^ ^(E)a(E) ^ f L\ f kL\ 
The same weight serves well in both the sine transform and also the nonlinear transforms of 
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Figure 3: Here we see the weighted Fourier transform of the spectrum without oscillations 
(blue solid curve) and with oscillations in the case of the normal (black solid curve) and 
inverted (red dashed curve) hierarchies. One can see that the solid, blue unoscillated curve 
is very close to zero. We have checked that this curve is essentially independent of the cutoff 
and so the reactor spectrum no longer affects RL. Comparing with Fig. 1 one can see that 
the difference RL between the depths of the minima is even greater in this weighted case, 
allowing for a better determination of the hierarchy than was possible with an unweighted 
Fourier transform. 

Ref. [10] which determine the hierarchy more reliably than RL+PV at baselines below about 
55 km [24]. As can be seen by comparing the unweighted and weighted cosine transforms in 
Figs. 1 and 3, not only does the weighting procedure preserve RL, but it actually increases 
the difference in the peak sizes between the normal and inverted hierarchies. Thus this 
solution to the dependence upon the high energy neutrino tail not only removes the spurious 
dependence, for any high energy reactor spectrum, but it increases the chance of success of 
the determination of the hierarchy. In simulations [24] we will show that when the weight 
function is optimized with a neural network this increase is of order 2%. 
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